Viability of Am^ ^ 1 eV^ sterile neutrino mixing models in light of MiniBooNE 
electron neutrino and antineutrino data from the Booster and NuMI beamlines 



G. KaragiorgiiQ Z. Djurcic^lll J. M. ConradiH M. H. Shaevitz^l and M. Sorel^l 

^Massachusetts Institute of Technology, Cambridge, MA 02139 
^Department of Physics, Columbia University, New York, NY 10027 and 
^Instituto de Fisica Corpuscular, IFIC, CSIC and Universidad de Valencia, Spain 



o 



in 

^\ 

in- 

m 
> 

On 

O 

On 



(Dated: May 7, 2010) 

This paper examines sterile neutrino oscillation models in light of recently published results from 
the MiniBooNE Experiment. The new MiniBooNE data include the updated neutrino appearance 
results, including the low energy region, and the first antineutrino appearance results, as well as 
first results from the off-axis NuMI beam observed in the MiniBooNE detector. These new global 
fits also include data from LSND, KARMEN, NOMAD, Bugey, CHOOZ, CCFR84, and CDHS. 
Constraints from atmospheric oscillation data have been imposed. We test the validity of the 
three-active plus one-sterile (3+1) and two-sterile (3+2) oscillation hypotheses, and we estimate 
the allowed range of fundamental neutrino oscillation parameters in each case. We assume CPT- 
invariance throughout. However, in the case of (3+2) oscillations, CP violation is allowed. With the 
addition of the new MiniBooNE data sets, there are clear incompatibilities between neutrino and 
antineutrino experiments under a (3+1) oscillation hypothesis. A better description of all short- 
baseline data over a (3+1) is provided by a (3+2) oscillation hypothesis with CP violation. However, 
we still find large incompatibilities among appearance and disappearance experiments, consistent 
with previous analyses, as well as incompatibilities between neutrino and antineutrino experiments. 
Aside from LSND, the data sets responsible for this tension are the MiniBooNE neutrino data set, 
CDHS, and the atmospheric constraints. On the other hand, fits to antineutrino-only data sets, 
including appearance and disappearance experiments, are found significantly more compatible, even 
within a (3+1) oscillation scenario. 

PACS numbers: 14.60.Pq, 14.60. St, 12.15.Ff 



I. INTRODUCTION 



Sterile neutrino oscillation models were proposed more 
than a decade ago as an explanation for the LSND 
anomaly [l|-[i|, an excess of events consistent with 
Ue oscillations at high Am^. These models relate Ue 
appearance (i^^ Ve) with and Ve disappearance 
(^/i ^// and ^^), motivating combined fits in all 

three oscillation channels. Relatively early in the discus- 
sion of models, it was demonstrated [5, 6] that a three- 
active plus one-sterile (3+1) neutrino oscillation model 
could not reconcile the LSND result with existing null 
results from other short-baseline (SBL) experiments, in- 
cluding KARMEN 0, NOMAD i, Bugey Q, CHOOZ 
[io|, CCFR84 ^ill, and CDHS [l2|, which had similar 
high Am^ sensitivity. However, it was shown that a 
three-active plus two-sterile neutrino (3+2) oscillation 
scenario provided a better description of these data sets 



In 2001, the MiniBooNE experiment began running 
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with the goal to test the LSND result using both neutrino 
and antineutrino beams. This is a short-baseline appear- 
ance and disappearance experiment located at Fermi Na- 
tional Accelerator Laboratory (Fermilab). MiniBooNE's 
first results, reported in 2007, described a search for 
^fi oscillations [l3|. These results were incom- 

patible with a simple two-neutrino oscillation interpre- 
tation of the LSND signal and, within this model, Mini- 
BooNE excluded the LSND result at the 98% CL. How- 
ever, this same analysis reported a 3.7cr excess of elec- 
tron neutrino candidate events at low energies, between 
300-475 MeV, which remains unexplained. Reference [l^ 
included the MiniBooNE first result in a global fit to all 
SBL experiments under the (3+1) and (3+2) oscillation 
scenarios. The analysis built on an earlier study, which 
introduced the possibility of CP violation (P(z^/x ^ ^e) 7^ 
^e)) within (3+2) fits p^]. Including the first 
MiniBooNE results into the global fit led to two obser- 
vations in Ref. Q: 1) MiniBooNE, LSND, and the null 
appearance experiments (KARMEN and NOMAD) are 
compatible under a (3+2) sterile neutrino oscillation sce- 
nario with large CP violation. 2) There is severe tension 
between appearance and disappearance experiments, at 
a level of more than 3a. In this paper we will consider 
both observations in light of new appearance data. Also, 
we will show that the incompatibility between appear- 
ance and disappearance experiments arises mainly from 
two disappearance data sets: CDHS and atmospheric 
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Data Set 


Channel 








LSND 


- 




BNB-MB(i/) 


- 




BNB-MB(f/) 


- 


Z^e 


NUMI-MB 


- 


Z^e 


KARMEN 


- 




NOMAD 


- 


Z^e 


Disappearance experiments: 






Bugey 


Z^e - 


■> 


CHOOZ 


Z^e - 


■> zy^ 


CCFR84 




■> zy^ 


CDHS 




■> z/^ 



TABLE I: Short-baseline oscillation data sets considered in 
this paper, and oscillation channel that each data set con- 
strains. 



constraints. 

Motivated by three new results from MiniBooNE, this 
paper re-examines the (3+1) and (3+2) global fits to 
the SBL data. These new results are: 1) an updated 
Yfi result 0; 2) first results for a ^ z^e search 

[17|; and 3) first ^/x results from the NuMI off- axis 

beam at MiniBooNE [18| . We consider these new results 
in combination with seven SBL data sets. These pro- 
vide constraints on: disappearance (from the CCFR84 
and CDHS experiments), disappearance (from the 
Bugey and CHOOZ experiments), Vy, z^e oscillations 
(from the NOMAD experiment), and v^^ oscilla- 

tions (from the LSND and KARMEN experiments). Fur- 
thermore, we have taken into account atmospheric con- 
straints based on the analysis of Ref. [l9|. These con- 
straints have been incorporated in our analysis following 
the method described in Ref. [l5[, and are included in 
fits to all SBL experiments, null SBL experiments, or as 
explicitly stated. Table |T] summarizes all SBL data sets 
used in the fits presented in this paper. 

In this work, we do not discuss experimental con- 
straints on sterile neutrino models other than SBL and 
atmospheric neutrino ones. Constraints from the mea- 
surement of the electron spectrum near the endpoint 
in beta-decay experiments are relatively weak as long 
as the mostly-sterile mass states are heavier than the 
mostly- active ones, because of the small electron flavor 
of the former (see Refs. [El, [20|). We make this assump- 
tion throughout the paper, by requiring that the heav- 
ier sterile neutrino mass eigenstates, and m4, obey 
ms > m4 > mi. Constraints on the energy density 
(and mass) in the Universe carried by sterile neutrinos 
from cosmic microwave background, matter power spec- 
trum, and supernovae data have been studied in Ref. [21| . 
While relevant, these constraints are found to be weaker 
than SBL ones, since sterile neutrino states do not nec- 
essarily feature thermal abundances at decoupling. Con- 
straints on the number of relativistic degrees of freedom 
from the observations of cosmological abundances of light 



elements produced at the epoch of Big Bang Nucleosyn- 
thesis may also be relevant, although model-dependent. 
For such a study involving one sterile neutrino species 
participating in the mixing, see for example Ref. [22|. 

The paper is organized as follows. In Section [III we 
provide a short description of the MiniBooNE experi- 
ment and the new data sets. In Section Hill we specify 
the formalism used in this analysis to describe (3+n) os- 
cillations, where n is the number of sterile neutrinos. In 
Section [IVl we discuss the analysis method followed, and 
describe in detail the way in which the three MiniBooNE 
data sets have been incorporated. In Section |Vl we 
present the results obtained for the (3+1) (CP-conserving 
only), and (3+2) CP-conserving and CP-violating hy- 
potheses. For each hypothesis, we quote the compati- 
bility between various sets of SBL experiments and re- 
port the best-fit neutrino mass and mixing parameters 
derived from the combined analysis of all experimental 
data sets. In the (3+2) CP-violating case, we discuss 
the constraints on the CP violation phase, inferred from 
a combined analysis of all SBL oscillation results. Fi- 
nally, in Section |Vll we discuss constraints to the (3+2) 
CP-violating models from each of the SBL experiments 
considered in this analysis. The goal of this particular 
study was to investigate whether the source of tension 
between appearance and disappearance experiments Q 
is a result of a single experiment, other than LSND. 



II. THE NEW MINIBOONE DATA SETS 

The MiniBooNE experiment uses a muon neutrino 
beam produced by 8 GeV protons from the Fermilab 
Booster Neutrino Beamline (BNB) impinging on a beryl- 
lium target. The target is located within a magnetic fo- 
cusing horn [23I . The current of the horn can be reversed 
for running neutrinos or antineutrinos, allowing Mini- 
BooNE to perform both neutrino and antineutrino oscil- 
lation searches. The detector is located L = 541 m 
from the primary target, and the neutrino flux has an 
average energy of ~ 0.75 GeV. This design maintains the 
LSND L/E oi ^ 1 m/MeV. The detector consists of a 
spherical tank with a 610-cm active radius, instrumented 
with 1520 8-inch photomultipliers. This is filled with 800 
tons of pure mineral oil. An outer veto region rejects cos- 
mic rays and neutrino events producing particles which 
cross the detector boundaries. 

The MiniBooNE neutrino data set used in this analysis 
corresponds to the updated results recently reported by 
the MiniBooNE collaboration [l6| . Compared to the first 
MiniBooNE result which was released in 2007 [III, the 
new result involves a re-analysis of the MiniBooNE low 
energy excess events and several updates to the Monte 
Carlo prediction. These updates include a new model 
of photonuclear effects, incorporation of new data on tt^ 
production and a better treatment of pion re-interaction 
in the detector following decay, an improved estimate and 
rejection method of the background from interactions 
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outside the detector, and improvements to the deter- 
mination of systematic errors. The updated low-energy 
analysis has resulted in a reduction to the significance 
of the excess from 3.7<j in the original analysis to 3.4cr, 
along with some slight modification to the shape of the 
energy spectrum; specifically, the peak of the excess has 
shifted slightly to higher neutrino energies. In addition, 
the new analysis extends in energy down to 200 MeV, 
compared to 300 MeV in the original analysis, which of- 
fers additional L/ E information. The new result also cor- 
responds to modestly higher statistics, corresponding to 
the total data collected during the experiment's neutrino 
running of 6.46 x 10^^ protons on target (POT), compared 
to 5.58x10^^, previously. 



More recently, the MiniBooNE Collaboration reported 
its first results from a search for oscillations, 
using a muon antineutrino beam [17]. The antineutrino 
analysis performed by MiniBooNE mirrors the updated 
neutrino analysis [16], and includes the Monte Carlo pre- 
diction updates of the latter. The total antineutrino data 
set used in the analysis corresponds to 3.39x10^^ POT. 
However, due to meson production and cross-section ef- 
fects, the antineutrino event sample, unlike the neutrino 
event sample, is statistically limited. Unlike the neutrino 
search, the MiniBooNE antineutrino search provides a 
direct test of the LSND result, similar to the search 
performed by KARMEN. The MiniBooNE sensitivity to 
^/i extends into the low-Am^ region allowed by a 

combined analysis of KARMEN and LSND data. Never- 
theless, the MiniBooNE antineutrino search has observed 
no conclusive signal, and a limit has been set, which is 
considerably weaker than the sensitivity, and comparable 
to the KARMEN limit. The limit degradation with re- 
spect to the sensitivity is due to a 2.8cr fluctuation of data 
above expected background observed in the 475-675 MeV 
energy region. Thus, at present, the MiniBooNE antineu- 
trino result is inconclusive with respect to oscillations al- 
lowed by LSND. However, MiniBooNE is in the process 
of collecting more antineutrino data. This is expected 
to improve the experiment's sensitivity to oscil- 
lations. Updated results are expected after about three 
years of running. 



The third new data set ^| arises from the fact that the 
MiniBooNE detector is illuminated by the off- axis (110 
mrad) neutrino flux from the NuMI beamline at Fermi- 
lab. This analysis has reported a 1.2cr excess of i^g-like 
events in the neutrino energy range below 900 MeV. The 
NuMI data set corresponds to a mean LjE that is ap- 
proximately the same as those of the MiniBooNE and 
LSND data sets, and therefore probes the same Am^ 
range, providing complementary information in oscilla- 
tion fits with MiniBooNE and LSND. 



III. (3+n) STERILE NEUTRINO OSCILLATION 
FORMALISM 

The formalism used in this paper follows that which 
was presented in Ref. [l5|. We provide a brief summary 
here. 

In sterile neutrino oscillation models, under the as- 
sumptions of CPT invariance and negligible matter ef- 
fects, the probability for a neutrino produced with flavor 
OL and energy to be detected as a neutrino of flavor /3 
after travehng a distance L, is given by [25, 26]: 

Piyo. ^ y^) = - 4 E.>, n{u:^^u^^Ua.jU;^) sin^ + 

2 E^>J im.U^^UajU;^) Sin 2x,j (1) 

where 1Z and 1 indicate the real and imaginary parts 
of the product of mixing matrix elements, respectively; 
a, /3 = e^ji^T^ or s, {s being the sterile flavor); i,j = 
1,...,3 + n (n being the number of sterile neutrino 
species); and xij = 1.27Amf^L/£^. In defining x^j, we 

take the neutrino mass splitting Amf^- = m'^ — m'j in eV^, 
the neutrino baseline L in km, and the neutrino energy 
E in GeV. For antineutrinos, the oscillation probability 
is obtained from Eq. [T] by replacing the mixing matrix 
U with its complex-conjugate matrix. Therefore, if the 
mixing matrix is not real, neutrino and antineutrino os- 
cillation probabilities can differ. 

For 3+n neutrino species, there are, in general, 2 + n 
independent mass splittings, (3 + n)(2 + n)/2 indepen- 
dent moduli of parameters in the unitary mixing matrix, 
and (2 + n)(l + n)/2 Dirac CP- violating phases that may 
be observed in oscillations. In SBL neutrino experiments 
that are sensitive only to i'^ ^ fyi^ I'e ^ and Ve 
transitions, the set of observable parameters is reduced 
considerably. In this case, the number of observable pa- 
rameters is restricted to n independent mass splittings, 
2n moduli of mixing matrix parameters, and n — 1 CP- 
violating phases. Therefore, for (3+2) sterile neutrino 
models (n = 2 case), for example, there are two indepen- 
dent mass splittings, Aml^ and /S.mli^ both defined to 
be greater than zero, four moduli of mixing matrix pa- 
rameters I UeA I , I ^/i4 1 , I \ , | U \ , and one CP-violating 
phase. The convention used for the CP-phase is: 

^45 = arg{U;,Ue,U^4U:^). (2) 

In that case, the general oscillation formula in Eq. [1] be- 
comes: 

Pil^a -^iya) = l- 4[(1 - |t/c.4|' - \Ua5\^) ■ 

(|/7c4|^sin^ X41 + |t/a5|^sin^ ^51) + 

\Ua4\^\Ua5\^Sm^X,^] (3) 

and 

PiUa y^^a) =4|[/c,4p|/7/34psin^ X41 + 

4|/7«5|^|t//35|^sin2 X51 + 

8|^a5||^/35||^a4||^/34| siu X41 siu X51 COs(x54 - (/)45) (4) 
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The formulas for antineutrino oscillations are obtained 
by substituting (j)^^ —045- 

For the case of (3+1) sterile neutrino models (n = 
1 case), the corresponding oscillation probabilities are 
obtained from Eqs. [3]and|4]by setting X51 = X54 = and 
\Ua5\ = 0. Note that, under the above assumptions, no 
CP violation is allowed for (3+1) models. 



IV. ANALYSIS METHOD 

In this section, we first provide an overview of the fit- 
ting technique. We then focus on the method followed 
for including the MiniBooNE data sets. The physics- 
and statistical-assumptions for the other null SBL exper- 
iments and LSND, which are also included in the fit, are 
described in detail in Ref. [5]. The constraints from at- 
mospheric experiments, according to Ref. [19], have been 
incorporated as described in Ref. [15l] . 



A. General Technique 

The Monte Carlo method used to apply the oscillation 
formalism described in Section HITl closelv follows the one 
described in Ref. [l5|. We start by randomly varying 
sets of oscillation parameters: Aml^^, |t/e4|, I^mI foi" the 
case of (3+1); Am|i, |/7e4|, |/7^4|, Am|i, |[/e5|, |f/^5|, (?^45 
for the case of (3+2). Without loss of generality, we 
take Amli > Amii. In CP-conserving models, ^45 is 
set to or TT by default, whereas in CP- violating models 
045 is allowed to vary within the full (0, 27r) range. For 
each set of oscillation parameters, a signal prediction is 
obtained and compared to observed data for each SBL 
experiment, in the form of a fo^ each experiment. For 
each set of oscillation parameters that is generated, the 
various x^'s are linearly summed together to form Xsbl^ 
which is then used to extract the best-fit parameters and 
allowed regions. 

A minimization is carried out using a Markov Chain 
p7| . This minimization procedure relies on calculating 
the difference between successive sets of parameters 
and using that as a measure of whether the new point in 
parameter space is a "good" point to step to, or whether 
a new set of parameters needs to be drawn again. This 
is realized in the form of a probability of accepting a new 
set of parameters, P{xi ^c^+i), given by 

P{xi Xi^i) = mm(l, e-(^^+i-^^)/^), (5) 

where Xi and x^+i are two successive points in param- 
eter space, and T is a Temperature parameter. Larger 
T values allow for larger Ax^ jumps on the x^ surface, 
and therefore by varying the T value, one can avoid lo- 
cal minima, as well as finely scan the parameter space. 
This minimization method is particularly preferred in fits 
with large parameter space dimensionality, as in the case 
of (3+2) oscillation fits, due to its higher efficiency. 



In extracting the various confidence level contours, we 
marginalize over the parameter space and report results 
obtained with Ax^ levels corresponding to 1 degree of 
freedom for exclusion limits, and 2 degrees of freedom 
for allowed regions. 

To quantify the statistical compatibility between vari- 
ous data sets, we use the Parameter Goodness-of-fit (PC) 
test introduced in [28]. In this test one quantifies how 
well various data sets are in agreement, by comparing 
the minimum x^ obtained by a fit where all data sets 
have been included as constraints, Xmin aii^ 
of the x^ minima obtained by independent fits for each 
experiment, i.e., 

(6) 

i 

where i runs over experimental data sets in consideration. 
The PC is obtained from x%g based on the number of 
common underlying fit parameters, ndfpG- 

PG=prob{xlG.ndfpG). (7) 

For example, for testing the compatibility between KAR- 
MEN and LSND for the (3+2) CP-conserving oscihation 
hypothesis, we fit for both KARMEN and LSND simul- 
taneously to extract Xmin k-\-l^^^^ KARMEN and 
LSND separately to extract xlrin,K^ and Xmin,L^ respec- 
tively, and obtain: 

)• (8) 

As these are appearance experiments, there are 4 com- 
mon fit parameters for a CP-conserving (3+2) model 
(Am|i, Am|i, \Ue4\\U^4\, and |t/e5| l^/xsl); therefore, 

PG{K, L) = probixlc {K,L)A)^ (9) 

It should be noted that x^-pi'obabilities and PC tests 
can lead to drastically different results [28| . This is often 
a consequence of a large data set simultaneously fitted 
with small data sets, where the large data set dominates 
the x^ of the simultaneous fit. 



B. Inclusion of the MiniBooNE Neutrino and 
Antineutrino Data Sets 

The MiniBooNE neutrino data set (BNB-MB(i/)), de- 
scribed in Secini is included in the fits in the form of two 
side-by-side distributions of Ve and charged-current 
quasi-elastic (CCQE) events. Each distribution is a func- 
tion of neutrino energy, reconstructed under the hypoth- 
esis of CCQE neutrino interaction kinematics, E^^ . The 
full 200-3000 MeV range of Ve CCQE data is used in the 
fit. The observed event distributions are compared to 
the corresponding Monte Carlo predicted distributions, 
and a x^ is calculated using a covariance matrix which 
includes systematic and statistical uncertainties as well 
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as systematic correlations between the predicted v^. and 
distributions. 

During the fit, we vary the v^^ distribution according 
to the sterile neutrino oscillation parameters, but keep 
the distribution unchanged. The Vy, distribution re- 
mains unchanged during the fit, despite the possibility of 
Vy, disappearance in the MiniBooNE data. In fact, Mini- 
BooNE has released results from and Vy, disappearance 
searches at Am^ 1 eV^ [29]. These results are relevant 
as constraints to sterile neutrino mixing parameters un- 
der consideration, but they have been purposefully omit- 
ted in this analysis, due to the fact that the and Vg 
CCQE samples used in the disappearance analysis [29| 
and the (different) and CCQE samples used as 
constraints in the appearance analyses [l6|, [Tg are highly 
correlated samples, and these correlations have not yet 
become available. We assume that including MiniBooNE 
Vy disappearance would have a small effect on sterile neu- 
trino fit results, given the large overlap of the disap- 
pearance limit from MiniBooNE with the corresponding 
limits from CDHS and CCFR84 [29^. However, the im- 
pact of the MiniBooNE disappearance results on the fits 
considered in this paper will be discussed. Nevertheless, 
we employ this side-by-side fitting method as it takes ad- 
vantage of correlations in the and z/g predictions and 
in order to effectively constrain the z^g prediction and re- 
duce systematic uncertainties in the z/g search. 

The fit method follows the details described in [l6| , ex- 
cept that it uses a different definition for the covariance 
matrix used in the calculation. Ref. involves an it- 
erative fit method where the calculation for each point 
on the parameter space being probed uses the covariance 
matrix calculated according to the best-fit signal predic- 
tion. Instead, in the MiniBooNE fits presented here, the 
surface is estimated using the covariance matrix cal- 
culated according to the signal prediction at each point 
of the parameter space under consideration. The two fit 
methods yield similar results, although, by definition, the 
iterative method of [l6| results in a relative shift of the 
allowed region to the left, i.e. towards smaller oscillation 
amplitudes. 

The MiniBooNE antineutrino data set (BNB-MB(z>)), 
described in Sec. [Ill is included in the fits in the same 
way as the BNB-MB(z^) data set, in the form of two side- 
by-side E^^ distributions of z>e and i>y CCQE events. 
In this case, the disappearance limit obtained using the 
MiniBooNE CCQE sample provides substantial cover- 
age of so-far unexplored sterile neutrino mass and mixing 
parameter space [29|. Even though we do not explicitly 
fit the MiniBooNE z/^ CCQE distribution for disappear- 
ance, we comment on the effect of the limit from [29| 
in Sec. |Vl and justify that excluding the MiniBooNE 
disappearance information from the fits does not sub- 
stantially affect the parameter space of interest. The full 
200-3000 MeV range of z/g CCQE data is used in the fit. 
The BNB-MB^^) data fit method also follows the details 
described in |17| , except that it uses the definition for the 
covariance matrix described above. 



In fits where both neutrino and antineutrino data are 
included, it has been assumed that the two data sets are 
fully uncorrelated. In reality, the two data sets have large 
systematic correlations. However, neglecting the effects 
of these correlations is a reasonable approximation, given 
that the antineutrino data set is statistics limited. 



C. Inclusion of the NuMI-beam Data Set 

The NUMI-MB data set [isj, described in Sec. M is 
included in the fits in the form of a distribution of z^e 
CCQE events as a function of reconstructed neutrino en- 
ergy, E^^ . The predicted z/g distribution is obtained by 
adding to the expected Ve CCQE background any con- 
tribution from Uy to z^e oscillations. The contribution 
is estimated as follows: First, a fully (100%) oscillated 
NUMI-MB Uy Ue sample is obtained by reweighting 
the BNB-MB fully oscillated Ue Monte Carlo pre- 

dicted sample according to the ratio of the NuMI-beam 
flux from [l^ to the BNB-MB flux [23|, as a function 
true neutrino energy. As the oscillation parameters vary 
during the fit, a signal prediction is calculated by rescal- 
ing the number of events in this fully oscillated sample 
by the corresponding oscillation probability, according 
to the true neutrino energy and distance travelled, from 
production to detection, of each event. We assume a 
constant L of 700 meters. The prediction is compared to 
the observed Ue CCQE events as a function of 10 bins of 
E^^ . The background and signal prediction are assumed 
to have the same fractional systematic uncertainties, and 
a statistical uncertainty is calculated for each point in 
the parameter space according to the signal prediction 
of each point under consideration. The data and back- 
ground central value and systematic uncertainty per 
bin have been estimated from [l8|. Unlike the system- 
atic uncertainties of the BNB-MB and z/g CCQE data 
sets, the NUMI-MB CCQE systematics have not been 
constrained using information from the CCQE spec- 
trum from the NuMI beamline. Furthermore, we have 
not considered potential systematic correlations among 
the z/g CCQE bins as a function of E^^ . 



V. (3-hl) AND (3+2) MODELS AFTER THE 
NEW MINIBOONE zy, AND NUMI RESULTS 

A. (3+1) FIT RESULTS 

In this section, the new MiniBooNE results are exam- 
ined under a (3+1) oscillation hypothesis and compared 
to LSND and other null SBL experiments. The new data 
sets are studied first within the context of appearance- 
only experiments, and subsequently in fits involving both 
appearance and disappearance data. Fits to only an- 
tineutrino and only neutrino SBL experiments are also 
explored. 
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TABLE II: Short-baseline oscillation fits considered in this paper. 



BNB-MB(v) BNB-MB(v) LSND 




-3 -2 -1 -3 -2 -1 -3 -2 -1 

10 10 10 10 10 10 10 10 10 



sin'(2e^e) sin2(2e^e) sin2(20^e) 

FIG. 1: Allowed regions (filled areas) at 90% and 99% CL from BNB-MB(zy)-only, BNB-MB(z/)-only, and LSND-only (3+1) 
fits. These fits are, by construction, CP-conserving. The stars indicate the three respective best-fit points. All three data sets 
show closed contours at 90% CL. See text for more details. 



Table iniprovid es a reference for all the different data 
set combinations explored in fits in this paper. 



1. Studies with appearance- only experiments 

Figure [1] shows the allowed regions obtained by in- 
dependent fits to each of the following three data sets: 
BNB-MB(i/), BNB-MB(i/), and LSND. The regions are 
estimated using a 2-dimensional global scan of the (3+1) 
parameter space (sin^ 2^^e7A^4i)- Each contour is 
drawn by applying a flat Ax^ = x^(sin^ 2^^e7 A^li) — 
Xmin ^^^^ surface, with respect to the global 

minimum returned by the fit. The figure shows that, 
similarly to LSND, both BNB-MB data sets yield con- 
tours which exclude the no-oscillations (null) hypothesis 
at 90% CL. The null x^'s correspond to 22.2 and 24.5 for 
BNB-MB (i/) and BNB-MB (i/), respectively. The closed 
contours reflect a contradiction to the oscillation results 
published by the MiniBooNE collaboration; this is a con- 
sequence of the different x^ definition involved in the fit 
method used here, as pointed out in Sec. IIVBI All three 



data sets, BNB-MB(i/), BNB-MB(i/), and LSND, yield 
similar best-fit parameters, indicated by the stars on the 
three graphs, of Am^^ of order a few eV^ and sin^ 2^^e 
of order 10~^ — 10~^. The minimum x^ and best-fit pa- 
rameters returned by each experiment are summarized in 
Table inD 

In light of the above BNB-MB results and the already 
established LSND anomaly, we find it interesting to study 
sterile neutrino oscillations with the LSND, BNB-MB (i^), 
and BNB-MB (i/) data sets assumed to be (positive) "sig- 
nal" experiments under both the (3+1) and the (3+2) 
models. This classification is based on the fact that all 
three data sets exclude the null result at 90% confidence 
level under a (3+1) oscillation hypothesis. Figure [2] shows 
the BNB-MB(i/) and BNB-MB(z/) event distributions for 
both the null and the best-fit (3+1) oscillation hypoth- 
esis for each data set. In the case of the BNB-MB (z/) 
data set, even though the best-fit hypothesis provides 
a better description of the event spectrum at 90% CL 
(Ax^ = xluii - Xlest-fit =4-7, for 2 fit parameters), it 
fails to fully explain the low energy excess. Therefore, 
the (3+1) oscillation hypothesis alone seems inadequate 
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Data Set 




X^-probability 


A 2 


sin^ 20^e 


sin^ 26^^ sin^ 20ee 


Appearance-only fits: 














LSND 


3.4 


(3) 


34% 


8.19 


0.0085 


- 


BNB-MB(zy) 


17.5 


(16) 


35% 


3.12 


0.0018 


- 


BNB-MB(zy) 


17.6 


(16) 


35% 


4.46 


0.0065 


- 


NUMI-MB 


2.0 


(8) 


98% 


6.97 


0.020 


- 


KARMEN 


6.0 


(7) 


54% 


6.81 


0.00096 


- 


NOMAD 


33.3 


(28) 


31% 


53.3 


0.00012 


- 


signal APP 


50.3 


(39) 


11% 


0.045 


0.98 


- 


signal APP* 


50.4 


(39) 


10% 


0.15 


0.090 


- 


null APP 


46.6 


(4-7) 


49% 


0.040 


1.00 


- 


APP 


97.1 


(88) 


24% 


0.045 


1.00 


- 


APP* 


97.2 


(88) 


24% 


0.15 


0.090 


- 


LSND + MB-BNB(i/) 


22.3 


(21) 


38% 


4.55 


0.0074 


- 


LSND + MB-BNB(iy)* 


22.3 


(21) 


38% 


4.55 


0.0074 


- 


LSND + MB-BNB(i/) + KARMEN 


33.6 


(30) 


29% 


0.57 


0.0097 


- 


Ri\rR A/rRrry'i -I- i\rTTA/rT a/tr -i- i\roA/rAD 


57.8 


(56) 


/O 


033 


1 00 




Appearance and disappearance fits: 














all SBL* 


197A 


(196) 


46% 


0.92 


0.0025 


0.13 0.073 


u 


90.5 


(90) 


47% 


0.19 


0.031 


0.031 0.034 


V 


87.9 (103) 


86% 


0.91 


0.0043 


0.35 0.043 



TABLE III: Comparison of best-fit values for mass splittings and mixing angles obtained from (3+1) fits to appearance data 
sets and appearance+disappearance data sets. Mass splittings are shown in eV^. The minimum from each fit, as well as 
the x^-probability are also given. The signal appearance (APP) data sets include BNB-MB(zy), BNB-MB(z/) and LSND. The 
nuh APP data sets include KARMEN, NOMAD and NUMLMB; the maximal best-fit sin^ 26>^e in this case is inconsequential, 
as it corresponds to a best-fit Am^ region of very poor sensitivity. See text for more details. 

* In these fits, the electron and muon content of the sterile neutrino mass eigenstate have been explicitly constrained to <0.3. 



2.5-, 



^ 1.5 



0.5 



T 



BNB-MB(v) 

data, statistical error 
null, systematic error 
(3+1) best-fit, systematic error 




^00 400 600 800 1000 1200 

E?^ (MeV) 



1500 



3000 



0.5 
0.45 

0.4 
0.35 

0.3 



i2 0.25 

I 0.2 
LU 

0.15 
0.1 
0.05 



BNB-MB(v) 

■ data, statistical error 
I null, systematic error 
I (3+1) best-fit, systematic error 



I I I I I 



J_ 



200 400 600 800 1000 1200 1500 
=QE, 



3000 



Ev"'= (MeV) 



FIG. 2: Left: Nuh and (3+1) best-fit predicted event distributions (Amli, sin^(26';,e)) = (3.12, 0.0018) for BNB-MB(zy). 
Right: Nuh and (3+1) best-fit predicted event distributions (Amli, sin^ (26*^6)) = (4.46, 0.0065) for BNB-MB(zy). The event 
distributions are shown as functions of reconstructed neutrino energy, E^^ . The data are shown in black points with statistical 
uncertainty. The null (no-oscillation) prediction is shown by the light gray histogram with (solid) systematic error band. The 
best-fit prediction (signal and background) is shown by the blue (dark gray) histogram with (shaded) systematic error band. 



as an explanation for the low energy excess, as also re- 
ported by the MiniBooNE collaboration [13, [TH. In the 
case of the BNB-MB(i^) data set, the best-fit hypothesis 
provides a better description of the data in the 500-700 
MeV range. However, the statistical uncertainties are too 
large to allow for a strong conclusion. 



The allowed regions obtained by a joint analysis of 
BNB-MB(i/) + BNB-MB(i/) + LSND, as well as a joint 
analysis of BNB-MB(i>) + LSND are shown on the left 
panels of Figs. [3] and IH respectively. In the case of the 
combined fit of all three data sets (Fig. [3]), due to the 
difference in preferred mixing amplitudes mentioned in 
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Data Set 


X^-probability (%) 


PG (%) 


APP 


24 


PG( BNB-MB(i/),BNB-MB(iy), LSND, NUMI-MB,KARMEN, NOMAD ) 


= Prob( 17.3,10 ) = 


6.8 


signal APP 


11 


PG( BNB-MB(zy),BNB-MB(i/),LSND ) 


= Prob( 11.9,4 ) = 


1.8 


LSND + MB-BNB(i/) 


38 


VCXi BNB-MBT/y'l T,SND 


— PrnbT 1 4 9 — 


49.0 


V APP 


29 


PG( BNB-MB(zy),LSND,KARMEN ) 


= Prob( 6.7,4 ) = 


15.3 


V APP 


40 


PG( BNB-MB(zy),NUMLMB, NOMAD ) 


= Prob( 4.9,4 ) = 


29.8 


all SBL* 


46 


PG( BNB-MB(zy),BNB-MB(iy),NUMLMB,LSND,KARMEN, 










NOMAD, Bugey,CHOOZ,CCFR84,CDHS,ATM ) 


= Prob( 42.0,18 ) = 


0.11 






PG( APP,DIS ) 


= Prob( 14.8,2 ) = 


0.06 






PG( v,v ) 


= Prob( 18.8,3 ) = 


0.03 


V 


47 


PG( BNB-MB(zy),NUMI-MB, NOMAD, CCFR84,CDHS, ATM ) 


= Prob( 14.7,8 ) = 


6.5 


V 


86 


PG( BNB-MB(zy),LSND,KARMEN,Bugey,CHOOZ ) 


= Prob( 8.43,7 ) = 


29.9 



TABLE IV: Summary of x^-probabilities for (3+1) fits with different combinations of SBL data sets, and PG results testing 
compatibility among different data sets. See text for more details. 

*/n these fits, the electron and muon content of the sterile neutrino mass eigenstate have been explicitly constrained to <0.3. 



the previous paragraph, the best-fit point ends up shift- 
ing from an intermediate Am^ and small mixing ampli- 
tude to a smaller Am^ and maximal mixing amplitude 
of 0.98. Obviously a maximal mixing amplitude is un- 
physical in the case of sterile neutrino oscillations. If the 
fits are repeated with the electron and muon content of 
the sterile mass eigenstate limited to values less than 0.3 
[38| . the returned x^-pi'obabilities are 10% and 38%, for 
BNB-MB(i/) + BNB-MB(i/) + LSND and BNB-MB(i/) 
+ LSND, respectively; the reduction in sin^ 2^^e space 
has essentially no effect on these results. The best-fit 
parameters from these fits are also given in Table IIIII 

Perhaps a more interesting observation regarding 
Fig.[T]is the striking similarity of BNB-MB(z/) and LSND 
90% CL allowed regions and best-fit oscillation param- 
eters, keeping in mind that both data sets describe an- 
tineutrino oscillations. It should be noted that in a (3+1) 
oscillation scenario, under the assumption of CPT invari- 
ance, there can be no difference between neutrino and an- 
tineutrino oscillation probabilities. However, a PG test, 
as described in Sec. IIV Al suggests a significantly higher 
compatibility (49%) between BNB-MB(i/) and LSND, 
rather than for all three signal experiments (BNB-MB(i^), 
BNB-MB(i/) and LSND) combined (1.8%). This is also 
supported by the x^-P^obabilities returned by the fits: 
11% in the case of the BNB-MB(z/) + BNB-MB(zy) + 
LSND fit, and 38% in the case of the BNB-MB(z/) + 
LSND fit. This incompatibility is due to the fact that 
the BNB-MB(z^) data set prefers a mixing amplitude ~3 
times smaller than the amplitude preferred by LSND or 
BNB-MB(i^), and excludes the LSND and BNB-MB(z/) 
best-fits at 99% CL. Table HVl provides a summary of the 
above x^-pi'obabilities and PG test results. 

Figures [3] and S] also illustrate the limits from various 
combinations of the remaining three (null) SBL appear- 
ance experiments, KARMEN, NOMAD, and NUMI-MB, 
under a (3+1) oscillation scenario, overlaid on the al- 
lowed regions described above. 

The 90% and 99% CL limits obtained by each of the 
null appearance experiments are shown on the left panel 



of Fig. [3l These limits correspond to the upper sin^ 20^e 
values allowed at each Am^^, estimated using a one- 
sided raster scan of the parameter space. It is interest- 
ing to point out that, despite the indication of a slight 
excess (1.2(t) of observed Ue-^ike events for neutrino en- 
ergies below 900 MeV found in the NuMI analysis [l8| , 
the currently assumed NUMI-MB systematic and statis- 
tical uncertainties are quite large, resulting in a limit that 
is much weaker relative to the limits of KARMEN and 
NOMAD. In fact, due to this excess and the large sys- 
tematic uncertainties, the NUMI-MB data set provides 
a very good fit to (3+1) models, with a x^-probability 
of 98%. The event distributions for the null and best- 
fit (3+1) oscillation hypothesis for the NUMI-MB data 
set are shown in Fig. [5l The observed distribution fits 
nicely to an oscillation signal at (Am^i, sin^(26>^e)) = 
(7.36, 0.019). Such large signal, however, would be in 
disagreement with the BNB-MB(z/) results. Additional 
data and reduced systematic uncertainties in the NUMI- 
MB analysis are necessary for higher sensitivity and more 
conclusive results. This is currently an ongoing effort 
and new results are expected soon. The limits from a 
combined NUMI-MB + KARMEN + NOMAD analysis 
are shown on the middle panel of Fig. [3l Both panels 
illustrate that the null appearance experiments provide 
essentially no constraints to the parameter space allowed 
by the BNB-MB and LSND data sets, except at higher 
Am^. 

The best-fit parameters obtained independently from 
the NUMI-MB and KARMEN data sets, shown in Ta- 
ble [ml are similar to those of LSND, BNB-MB(z/), and 
BNB-MB(i^). The NOMAD data set, on the other hand, 
prefers a much larger Amli ~50eV^, and a much smaller 
sin^ 20^e ^10-^. 

A combined analysis of all appearance data yields a 
X^-probability of 24% for the best-fit hypothesis, both in 
the case where maximal mixing is allowed in the fit, and 
in the case where the electron and muon content of the 
sterile mass eigenstate has been limited to small values 
(<0.3). The allowed region obtained by a joint analysis 
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FIG. 3: Left: Allowed 90% and 99% CL regions (light and dark filled areas, respectively) from a combined analysis of BNB- 
MB(z/), BNB-MB(zy) and LSND data sets, and 90% and 99% exclusion limits (light and dark curves, respectively) from 
each of the null appearance experiments, NUMI-MB (solid curves), KARMEN (dashed curves) and NOMAD (dotted curves). 
Middle: The same allowed region with overlayed 90% and 99% exclusion limits from a combined analysis of all null appearance 
experiments. Right: Allowed region obtained by a combined analysis of all appearance data sets, signal and null. See text for 
more details. 




FIG. 4: Left: The allowed 90% and 99% CL regions (light and dark filled areas, respectively) from a combined analysis of 
BNB-MB(z/) and LSND data sets, and 90% and 99% exclusion hmits (light and dark curves, respectively) from KARMEN. A 
comparison of only these three experiments is interesting, as these three experiments have searched for antineutrino oscillations 
at short baselines. Right: The allowed regions obtained from a combined analysis of all three experiments (BNB-MB(z>), LSND, 
and KARMEN). See text for more details. 



of all appearance experiments under a (3+1) oscillation 
scenario is shown in the right panel of Fig. [3l 

Similarly, Fig. |4] (left) corresponds to the allowed re- 
gion obtained by a joint analysis of BNB-MB(i/) + 
LSND. The limit shown is that of the KARMEN experi- 
ment, which is the only other SBL experiment to perform 
an appearance search with antineutrinos. The KARMEN 
limit provides substantial coverage of the joint LSND and 



BNB-MB(i^) allowed region, excluding the best-fit point 
of the LSND + BNB-MB(zy) fit at >99% CL. However, 
KARMEN imposes little constraint to the lower-Am^ al- 
lowed solutions. A joint analysis of all three data sets 
yields a x^-pi'obability of 29% for the best-fit hypoth- 
esis, and an allowed region shown on the right panel 
of Fig. m The x^-P^obability remains the same for fits 
where the electron and muon content of the sterile mass 
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FIG. 5: Null and (3+1) best-fit predicted event distributions 
(Amli,sin2(26';,e)) = (7.36, 0.019) for NUMI-MB. The data 
are shown in black points with statistical uncertainty. The 
null (no-oscillation) prediction is shown by the light gray his- 
togram with (solid) systematic error band. The best-fit pre- 
diction (signal and background) is shown by the blue (dark 
gray) histogram with (shaded) systematic error band. 



FIG. 6: The allowed 90% and 99% CL regions (light and dark 
filled areas, respectively) from a combined analysis of BNB- 
MB(z/) + BNB-MB(zy) + LSND data sets, and 90% and 99% 
exclusion limits (light and dark curves, respectively) from a 
combined analysis of all remaining (null, appearance and dis- 
appearance) SBL experiments. The null fit includes atmo- 
spheric constraints. The null SBL experiments exclude the 
joint 99% CL allowed region at 99% CL. 



eigenstates have been limited to values less than 0.3. As 
shown in Table |lVl the three data sets are compatible 
at 15.3%. New results from MiniBooNE with increased 
antineutrino statistics should be able to provide more in- 
formation to these fits [I^ • 



2. Studies with appearance and disappearance experiments 

Much stronger constraints than those of the null ap- 
pearance experiments are provided by the null disap- 
pearance experiments (CCFR84, CDHS, CHOOZ, and 
Bugey) and atmospheric constraints, under the assump- 
tions of CPT conservation and unitarity of the neutrino 
mixing matrix. The 90% and 99% CL exlusion limits 
from a combined analysis of all null data sets (NUMI-MB, 
KARMEN, NOMAD, Bugey, CHOOZ, CCFR84, CDHS, 
and atmospheric constraints) are shown in Fig. [6l The 
figure shows that the parameter space jointly- allowed by 
BNB-MB(i/) + BNB-MB(zy) + and LSND at 99% CL is 
excluded by a combined analysis of all null SBL experi- 
ments, appearance and disappearance, at 99% CL. The 
severe tension between LSND and the null SBL experi- 
ments [l3| continues to exist and in fact increases further 
with the addition of BNB-MB(i/) and BNB-MB(i^) data. 
The signal results show low (0.15%) compatibility with 
null results. The LSND result remains to be mostly re- 
sponsible for the low compatibility, as the BNB-MB(i^) 
and BNB-MB(i^) experiments show 14% and 3.7% com- 
patibility with the null experiments, respectively. 

Figure shows the allowed region obtained by the 
joint BNB-MB(i/) + LSND + KARMEN + Bugey + 
CHOOZ analysis. Here, the Ue disappearance constraints 



from Bugey and CHOOZ are interesting to consider from 
the perspective of a joint analysis of only antineutrino 
SBL experiments. In a joint fit, all of the above (an- 
tineutrino) experiments yield a high x^-P^'obability of 
86%, and 29.9% compatibility. In these fits, Bugey 
and CHOOZ constrain |t/e4|, but provide no direct con- 
straints on |t//x4|- However, a joint analysis with the 
LSND, BNB-MB(i^), and KARMEN appearance experi- 
ments, which are sensitive to the product of | /7e4 1 1 ^/24 1 , 
provides indirect constraints to the |/7^4| mixing element. 
Figure [3 (left) also shows that a fit to all antineutrino ex- 
periments without LSND yields similar closed contours 
at 90% CL, which include the best-fit point. Current 
constraints from MiniBooNE on disappearance [2^ 
provide relatively small constraints to the sin^ 26>^^ al- 
lowed space, as illustrated in the right panel of Fig. [71 
However, new results from a joint MiniBooNE and Sci- 
BooNE [30i] disappearance search, which are expected 
soon [3i[, may be able to probe this region with higher 
sensitivity, and will be interesting within the context of 
CPT-violating models. According to the best-fit oscil- 
lation parameters from a fit to only antineutrino SBL 
data, MiniBooNE should observe muon antineutrino dis- 
appearance with an amplitude of sin^ 2^^^ ~ 0.35, at 
Am^ - 0.91 eV^ The MINOS experiment [H should 
also have sensitivity to these oscillation parameters in an- 
tineutrino running mode; muon antineutrino disappear- 
ance search results from MINOS are expected soon [s^. 
Incorporation of the upcoming MiniBooNE and MINOS 
disappearance results in these fits is currently being in- 
vestigated. 

Neutrino-only fits also yield a reasonably high x^- 



11 



10" 



10e 



10-^ 



(3+1) 




V SBL 90% CL 
■v SBL 99% CL 

V SBL 3a CL 

□v SBL without LSND 90% CL 



10 



10 

s\n{2QJ 



10 



10" 



10 



> 

3. 1 

CM 5 

E 

< 

AO-' 



10"^ 



10 



(3+1) 



[vSBL 90% CL 
|v SBL 99% CL 

_MiniBooNE v„ 90% CL 



10 

sin'(2e^^) 



10 



FIG. 7: The allowed 90%, 99%, and 3a CL regions from a combined analysis of all antineutrino SBL data sets. The left plot also 
shows the 90% CL allowed region obtained from a combined analysis of all antineutrino experiments except LSND (KARMEN, 
BNB-MB(z/), Bugey, and CHOOZ). The right plot also shows the 90% CL exclusion limit from [29]. The MiniBooNE 
disappearance search excludes the parameter space to the right of the line at 90% CL. See text for more details. 




FIG. 8: The allowed 90%, 99%, and 3a CL regions from a combined analysis of all neutrino SBL data sets. The right plot also 
shows the 90% CL exclusion limit from The MiniBooNE disappearance search excludes the parameter space to the 

right of the line at 90% CL. See text for more details. 



probability of 47%; the corresponding allowed regions 
are shown in Fig. [51 Current constraints from Mini- 
BooNE disappearance are shown on the right panel 
of Fig. m Interestingly, fits to only neutrino SBL data 
also yield a closed contour at 90% CL. The parame- 
ter space, however, points to smaller mixing amplitudes 
relative to those preferred by the antineutrino-only fit. 
Neutrino-only fits and antineutrino-only fits are incom- 
patible, with a PC of 0.03%, as shown in Table |lVl The 
large incompatibility between antineutrino and neutrino 



SBL results suggests that the neutrino and antineutrino 
data sets cannot be accommodated within a (3+1) CPT- 
conserving sterile neutrino oscillation scenario. However, 
the constraining power of antineutrino SBL experiments 
alone on Am^^ and sin^ 2^^e is remarkable, and invites 
exploration of models that provide the possibility of dif- 
ferent oscillation patterns for neutrinos versus antineu- 
trinos. 

Figure [9] (left) shows a comparison of the BNB-MB(zy), 
BNB-MB(i^), and NUMTMB event distributions for the 
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FIG. 9: Left: MiniBooNE predicted event distributions using the neutrino-only (3+1) best-fit parameters 
(Am^, sin^(2^)) = (0.19, 0.031) in blue (dark gray) solid line and antineutrino-only (3+1) best-fit parameters 
(Am^, sin^(2^)) = (0.91, 0.0043) in blue (dark gray) dashed line. The null predictions are shown in light gray with system- 
atic error bands. The observed data are shown in black points with statistical error bars. Right: MiniBooNE predicted event 
distributions using the best-fit parameters obtained from a (3+2) CP- violating fit to all SBL experiments and appearance-only 
SBL experiments, in red (dark gray) solid line and red (dark gray) dashed line, respectively. 



neutrino-only best-fit parameters and antineutrino-only 
best-fit parameters. The neutrino best-fit parameters 
provide a better description to BNB-MB(i/) and NUMI- 
MB distributions than the antineutrino best-fit parame- 
ters, with xlNB-MBiu) = 18.4 and X%umi-mb = 4-4, 



compared to x]3nb-mb{v) = 32.4 and X%umi-mb = 
4.8. On the other hand, the antineutrino best-fit param- 
eters provide a better description to BNB-MB(i/) than 
the neutrino best-fit parameters {x'bnb-mb{p) ~ 19-7, 
compared to x\nb-mb{v) = 21.7). 
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FIG. 10: Allowed regions in (Amli,Amii) space for fits to CP-conserving (CPC, left) and CP- violating (CPV, right) (3+2) 
oscillation models. Only the BNB-MB(zy), BNB-MB(z/) and LSND data sets have been included in the fit. 



The best-fit results from the (3+1) oscillation fit in- 
volving all SBL data sets are summarized in Table 
lllli The best-fit parameters from neutrino-only and 
antineutrino-only fits are also shown. 



electron and muon content of the sterile mass eigenstates 
have been limited to values less than 0.3. This is a real- 
istic assumption for sterile neutrino oscillation models. 



B. (3+2) FIT RESULTS 

Neutrino oscillation models with more than one sterile 
neutrino have been of particular interest because they 
open up the possibility of observable CP violation effects 
in short-baseline neutrino oscillations. If (3+n) sterile 
neutrino oscillations are realized in nature, with n >1, 
CP violation becomes a natural possibility, which is very 
appealing from the perspective of theories attempting to 
explain the matter-antimatter asymmetry in our universe 

m 

In this section, the new MiniBooNE results are ex- 
amined under both a CP-conserving (CPC) and a CP- 
violating (CPV) (3+2) oscillation hypothesis. The new 
results are studied first within the context of appearance- 
only experiments, and subsequently in fits involving both 
appearance and disappearance data. 

From the point of view of the data at hand from LSND, 
BNB-MB(zy), and BNB-MB(zy) (see Fig. [I]), CP viola- 
tion offers the potential of reconciling two experimen- 
tal signatures — an excess in LSND data at 3.8a signif- 
icance and one suggested at 90% CL in BNB-MB an- 
tineutrino data, both pointing to relatively large mixing, 
reconciled with a possible excess found in BNB-MB neu- 
trino data suggesting relatively small mixing, both at a 
similar L/E — as manifestations of the same underlying 
oscillation hypothesis. 

It should be noted that in the studies presented in this 
section, due to the larger dimensionality of the fits, the 



1. Studies with appearance- only experiments 

Allowing for CP violation in (3+2) fits to LSND and 
BNB-MB (z/ and u) data leads to a significant reduction 
in absolute of 12.2, for 1 degree of freedom (dof), 
corresponding to a best-fit CPV phase ^45 =l.l7r. The 
X^-probability of the fit increases from 13% in the CPC 
case to 53% in the CPV case. The same test can be per- 
formed using all appearance data. In this case, allowing 
for CP violation leads to a reduction in of 13.3 for 
1 dof, with a best-fit CPV phase ^45 =l.l7r. The x^- 
probability from the CPV fit is comparable to that of a 
signal-only fit, at 56%. 

The 90% and 99% CL allowed (Am|i, Am|i) param- 
eter space obtained by a combined fit to BNB-MB (z/) 
+ BNB-MB (i/) + LSND is shown in Fig. [JOl The fig- 
ure illustrates that a CPV scenario (right panel) is much 
more restrictive in Am^ parameters compared to a CPC 
scenario (left panel). That is true both at 90% and 99% 
CL, shown by the significant reduction in allowed regions 
around Amli =0.5 eV^ and Amli =1 eV^. 

A similar effect is seen in the case of fits to all appear- 
ance experiments, as shown in Fig. [TTl Allowing for CP 
violation in fits to neutrino and antineutrino appearance 
data sets lead to a considerable improvement in the fit 
quality, and provides strong constraints to the Am^^ and 
Amli parameters of the model. 

The best-fit parameters for the signal-only and 
appearance-only fits are summarized in Table |Vl 
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FIG. 11: Allowed regions in (Amli,Amii) space for fits to CP-conserving (CPC, left) and CP-violating (CPV, right) (3+2) 
oscillation models. Only appearance data sets have been included in the fit. 
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fit. 



2. Studies with appearance and disappearance experiments 



A dramatic reduction in the allowed {Amli^ Amli) 
parameter space occurs once all SBL data sets are con- 
sidered in the fit, as shown in Fig. [121 Compared to 
the CPC hypothesis, with the addition of disappearance 
constraints, the CPV hypothesis fails to provide a sub- 
stantially better description of the data, reflected by the 
reduction in of Xcpc ~Xcpv =2.2 for 1 dof. Further- 
more, mainly due to CDHS [5], the allowed Am^ regions 
shift to higher Aml-^ values near Amli =25 eV^. 



The returned x^-probabilities from fits to all SBL data 
are 52% for the CPC fit, and 54% for the CPV fit. A PG 
test among all experimental data sets for the CPV case 
yields a compatibility of 7.0%. While the x^-probability 
and compatibility for the (3+2) CPV scenario are per- 
fectly acceptable, as will be discussed in Sec. |Vll an un- 
derlying source of tension exists due LSND and three 
other data sets: BNB-MB(i^), CDHS, and atmospheric 
constraints. The best-fit parameters extracted from a fit 
to all SBL data are also summarized in Table IVl 

A comparison of Tables Hill and IVl suggests that, with 
the addition of the new data sets from MiniBooNE, the 
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Data Set 


Fit 


(dof) 


X^-probability 


Ami, 




|t/e4| 


\U^4\ 


|t/e5| 


\U^5\ 


045 


signal APP 


CPV 


34.7(36) 


53% 


0.59 


1.21 


0.19 


0.33 


0.20 


0.16 


l.lTT 


signal APP 


CPC 


46.9(37) 


13% 


2.01 


2.22 


0.42 


0.24 


0.33 


0.33 





APP 


CPV 


82.5(85) 


56% 


0.39 


1.10 


0.40 


0.20 


0.21 


0.14 


l.lTT 


APP 


CPC 


95.8(86) 


22% 


0.18 


2.31 


0.32 


0.38 


0.086 


0.071 





all SBL 


CPV 


189.3(192) 


54% 


0.92 


26.5 


0.13 


0.13 


0.078 


0.15 


1.77r 


all SBL 


CPC 


191.5(193) 


52% 


0.92 


24.0 


0.12 


0.14 


0.070 


0.14 








TABLE V: Comparison of best-fit values for mass splittings and mixing parameters for (3+2) CP-conserving (CPC) and CP- 
violating (CPV) models. Mass splittings are shown in eV^. The appearance experiments include BNB-MB(zy and u), LSND, 
NUMI-MB, KARMEN, and NOMAD. The signal experiments include LSND, BNB-MB(zy), and BNB-MB(zy). See text for 
more details. 
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FIG. 13: Projection of Ax^ = ~ Xmin as a function of the 
CP- violating phase 045. The dashed horizontal lines indicate 
the 90% and 99% CL Ax^. 



(3+2) CPV oscillation hypothesis provides a better de- 
scription of all SBL data, compared to the (3+1) hypoth- 
esis. Compared to (3+1) models, (3+2) CP-conserving 
models give a reduction of 5.9 units for 3 additional 
fit parameters, while (3+2) CP-violating models give a 
reduction of 8.1 units with 4 additional parameters. 
This represents a relative improvement that is signifi- 
cantly smaller than that found in Ref. [5|] from fits using 
data sets prior to atmospheric constraints and the new 
MiniBooNE results. 

The MiniBooNE event distributions for the (3+2) CPV 
best-fit parameters are shown on the right panels of 
Fig. [9l The resulting x^-pi'obabilities are 8.6%, 6.7%, 
and 33.6%, for BNB-MB(i/), BNB-MB(i/), and NUMI- 
MB, respectively, obtained using the best-fit parameters 
from a (3+2) CPV fit to all SBL data. Notice, how- 
ever, that the best-fit parameters from a (3+2) CPV fit 
to appearance-only SBL data provide a better descrip- 
tion of all three MiniBooNE data sets than the best-fit 
parameters from a (3+2) CPV fit to all SBL data, par- 
ticularly for the BNB-MB(z/) data set. Furthermore, in 
the case of NUMI-MB, the (3+2) CPV appearance-only 



best-fit distribution, shown in dashed red (dark gray) on 
the right (x^-probability=61%), has comparable agree- 
ment with data as the (3+1) neutrino-only best-fit dis- 
tribution, shown in solid blue (dark gray) on the left 
(x^-probability=82%). In the case of BNB-MB(i/), the 
(3+2) CPV appearance-only best-fit parameters are ac- 
tually more preferred than the (3+1) neutrino-only best- 
fit parameters, with a x^-pi"ot)ability of 56.9%, rather 
than 30.1%. However, in the case of BNB-MB(i/) the x^- 
probability is highest (23.4%) for the (3+1) antineutrino- 
only best fit parameters. 

Figure [13] shows the projection of Ax^ = ~ Xmin 
as a function of the CP-violating phase ^45 for the 
three fits discussed in this section: the appearance-only 
fit projection is shown in the solid orange (light gray) 
line, the BNB-MB(i/)+BNB-MB(i/)+LSND fit projec- 
tion in dashed orange (light gray), and the projection 
from a fit to all SBL experiments is shown in blue (dark 
gray). Both the fit to the three signal experiments (BNB- 
MB(z/)+BNB-MB(i/)+LSND) and the fit to appearance- 
only experiments seem to strongly prefer a CPV phase at 
045 =l.l7r, as illustrated by the three overlapping dips in 
the Ax^ distribution. However, when fits to all SBL data 
are considered, the strong dependence disappears and a 
CPV phase at ^45 =1.77r is preferred. 



VI. CONSTRAINTS TO (3+2) CP- VIOLATING 
FITS FROM EACH SBL EXPERIMENT 

In this section we study the constraints to experimen- 
tally allowed (3+2) CP-violating oscillations by each of 
the SBL experiments considered in our fits. This is ac- 
complished through a study where fits are performed 
using all-but-one experiment at a time. Within this 
study, we are also interested in examining the source of 
incompatibility between appearance and disappearance 
data, as well as testing compatibility between neutrino 
and antineutrino appearance search results within a CP- 
violating scenario. The latter is motivated by the large 
incompatibility found in neutrino versus antineutrino fits, 
as well as appearance versus disappearance fits. 

Table IVII summarizes the x^-pi"ot)ability and PG re- 
sults from (3+2) CP-violating fits. The upper rows 
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Data Set 




X^-probability (%) 


PG (%) 


all SBL 


189.3 (192) 


54.2 


PG( BNB-MB(iy),BNB-MB(iy), NUMI-MB, LSND, KARMEN, 














NUMAD, Bugey, CHUUZ,UL/1^ R84,UDHS,A1 M ) 


= Prob( 53.9,40 ) = 


7.0 








v>c^ ( 




= Prob( 25.5,4 ) = 


0.004 








T>r^ ( 




= Prob( 25.4,7 ) = 


0.06 


APP 


82.5 (85) 


55.7 


PG( 


BNB-MB(zyj,BNB-MB(z^ j,N UMl-MB, LSND, KARMEN , 
















= Prob( 20.2,25 ) = 


73.6 


DIS 


81.3 (103) 


94.4 


v>o ( 


Jjugey, v_^rli^wZj,v^v_^Jr JrLo4t, v_yiJrlo, A j 


= Prob( 8.14,11 ) = 


70.1 


V 


81.3 (86) 


62.4 




tDi>i iD-NLiDl^h' j,iN U iVii-iViiD,iN \^iViAiJ,L7L7r xlo4,L7iJrlo, Ai iVi ) 


= Prob( 17.3,17 ) = 


43.4 


V 


82.6 (99) 


88.3 




iDi\ iD-NLiDl^h' j ,ivAxliVirL/iN ,ijoiN iJ,iDUgey,L7JrH^\^Zi ) 


= Prob( 11.2,16 ) = 


79.7 


V APP 


45.1 (53) 


77.0 


PG( 


BNB-MB(i/j,N UM1-MB,N(JMAD ) 


= Prob( 3.07,10 ) = 


98.0 


V APP 


27.1 (27) 


46.0 


PG( 


JDiNlD-MlD(^'j,KAKMEiN ,L/biNJJ ) 


= Prob( 6.88,10 ) = 


73.7 








PG( 


u AFF,!/ AFF ) 


= Prob( 10.3,5 ) = 


6.8 


all - BNB-MB(z/) 


167.3 (174) 


62.8 


PG( 


all SBL - BNB-MB(i/) , BNB-MB(i/) ) 


= Prob( 15.7,5 ) = 


0.78 


all - BNB-MB(i/) 


167.4 (174) 


62.6 


PG( 


all SBL - BNB-MB(i/) , BNB-MB(i^) ) 


— PrnbT 8 69 5 — 


13 


all - NUMI-MB 


183.7 (182) 


45.1 


PG( 


all SBL - NUMLMB , NUMLMB ) 


— PrnbT 3 90 5 — 


56 


all - LSND 


175.2 (187) 


72.2 


PG( 


all SBL - LSND , LSND ) 


— PrnbT T 9 5 5 — 


2.9 


all - KARMEN 


179.4 (183) 


56.1 


PG( 


all SBL - KARMEN , KARMEN ) 


= Prob( 4.53,5 ) = 


48 


all - NOMAD 


153.2 (162) 


67.8 


PG( 


all SBL - NOMAD , NOMAD ) 


= Prob( 1.96,5 ) = 


86 


all - Bugey 


140.4 (132) 


29.2 


PG( 


all SBL - Bugey , Bugey ) 


= Prob( 3.90,4 ) = 


42 


all - CHOOZ 


179.9 (178) 


44.6 


PG( 


all SBL - CHOOZ , CHOOZ ) 


= Prob( 3.09,4 ) = 


54 


all - CCFR84 


174.3 (174) 


47.9 


PG( 


all SBL - CCFR84 , CCFR84 ) 


= Prob( 0.35,4 ) = 


99 


all - CDHS 


172.8 (177) 


57.5 


PG( 


all SBL- CDHS , CDHS ) 


= Prob( 9.21,4 ) = 


5.6 


all - ATM 


184.0 (190) 


60.9 


PG( 


all SBL - ATM , ATM ) 


= Prob( 5.31,1 ) = 


2.1 



TABLE VI: Comparison of x^-probabilities for (3+2) CP- violating fits with different combinations of SBL data sets. Also shown 
are PG results testing compatibility among different data sets. The last eleven rows of the table provide the compatibility (PC) 
between the experiment being removed from each fit and all remaining experiments. See text for more details. 



Data Sets 


PG (%) 


Data Sets 


PG (%) 


APP vs. DIS 




0.004 


u vs. u 


0.06 


APP (no BNB-MB(i/)) vs. DIS (no 


CDHS + ATM) 


23.7 


V (no BNB-MB(i/) + CDHS + ATM) vs. v 


56.5 


APP (no BNB-MB(ry)) vs. DIS (no 


CDHS) 


0.36 


V (no BNB-MB(ry) + CDHS) vs. v 


3.7 


APP (no BNB-MB(i/)) vs. DIS (no 


ATM) 


0.52 


V (no BNB-MB(iy) + ATM) vs. v 


4.4 


APP (no BNB-MB(i/)) vs. DIS 




0.067 


V (no BNB-MB(zy)) vs. v 


1.1 


APP vs. DIS (no CDHS + ATM) 




2.9 


V (no CDHS + ATM) vs. v 


2.3 


APP vs. DIS (no CDHS) 




0.027 


V (no CDHS) vs. v 


0.07 


APP vs. DIS (no ATM) 




0.019 


V (no ATM) vs. v 


0.21 



TABLE VIL Comparison of compatibility between appear- 
ance (APP) and disappearance (DIS) experiments, within 
a (3+2) CP-violating scenario. The BNB-MB(zy) data set, 
CDHS data set, and atmospheric constraints (ATM) are re- 
moved from the fits as specified in order to establish the 
source of tension between appearance and disappearance ex- 
periments. Compatibilities are obtained using ndfpo =4. See 
text for more details. 



TABLE VIII: Comparison of compatibility between neutrino 
(ly) and antineutrino (p) experiments, within a (3+2) CP- 
violating scenario. The BNB-MB(zy) data set, CDHS data set, 
and atmospheric constraints (ATM) are removed from the fits 
as specified in order to establish the source of tension between 
neutrino and antineutrino experiments. Compatibilities are 
obtained using ndfpc =7. See text for more details. 



summarize x^-P^'obabilities and PCs from fits to all 
SBL experiments, as well as fits to appearance-only, 
disappearance-only, neutrino-only, antineutrino-only, 
neutrino appearance-only, and antineutrino appearance- 
only data sets. Appearance and disappearance data sets, 
as well as neutrino and antineutrino data sets, are in- 
compatible with a PG of less than 0.1%. Grouping SBL 
appearance-only data sets according to whether they are 
neutrino or antineutrino experiments yields significantly 
higher compatibilities — 98% among u appearance exper- 



iments, and 74% among u appearance experiments. The 
compatibility between u and u appearance-only results is 
lower, at 6.8% but still acceptable. In the case where dis- 
appearance experiments are included in the comparison 
between neutrino and antineutrino fits, the compatibil- 
ity among all u SBL data sets remains considerably high, 
at 80%, as does the compatibility among all u SBL data 
sets, at 43%. However the compatibility between u and 
u results is only 0.06%. 

The remaining rows of Table IVII provide the x^- 
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probabilities of global fits under the same oscillation sce- 
nario where one experiment is excluded from the fit at a 
time (as indicated by the sign in the table). The 
probabilities of all fits are acceptable, ranging between 
29.2% for a fit excluding the Bugey data set, and 72.2% 
for a fit excluding the LSND data set. Aside from LSND, 
three experiments stand out as having the poorest com- 
patibility when compared to a global fit with all other 
SBL data sets: 1) BNB-MB(zy), 2) CDHS and 3) atmo- 
spheric constraints (ATM). These three experiments have 
been identified as the possible source of tension between 
appearance and disappearance experiments, or neutrino 
and antineutrino experiments. The remaining combi- 
nations yield reasonably high compatibilities of at least 
42%, with the exception of LSND and BNB-MB(i/) which 
are compatible with the remaining data sets at 2.9% and 
13%, respectively. 

To further test the hypothesis that the tension between 
appearance and disappearance experiments is a result of 
the BNB-MB(i/) and CDHS data sets and atmospheric 
constraints, the compatibility between appearance and 
disappearance experiments is re-evaluated several times. 
Each time, a different combination of these three exper- 
iments is excluded from the fits. The results are sum- 
marized in Table IVIII The compatibility among appear- 
ance and disappearance experiments with BNB-MB(z/), 
CDHS, and atmospheric constraints excluded from the 
fits is high, at 23.7%. The BNB-MB(z/) data set alone 
is not responsible for the disagreement between appear- 
ance and disappearance experiments, as suggested by the 
sixth row of the Table IVH, Even with BNB-MB(z/) in- 
cluded in the fit, a compatibility of 2.9% can be obtained 
if CDHS and atmospheric constraints are excluded from 
the fit. 

The same test can be performed between neutrino and 
antineutrino experiments. The results are summarized 
in Table IVIIII Again, the compatibility between neu- 
trino and antineutrino experiments is re-evaluated sev- 
eral times; each time, a different combination of the 
BNB-MB(z/), CDHS, and atmospheric constraint data 
sets is excluded from the fits. Here, the compatibil- 
ity among neutrino and antineutrino experiments with 
BNB-MB(i/), CDHS, and atmospheric constraints ex- 
cluded from the fits is even higher, at 56%. The BNB- 
MB(i/) data set is just as responsible for the disagree- 
ment between neutrino and antineutrino experiments as 
the CDHS data set and atmospheric constraints alone. 
The tension seems to be caused by all three experiments, 
as none of them independently excluded from the fit can 
accound for the increase in compatibility from '^1% (or 
less) to 56%. 

It is possible that higher compatibility between BNB- 
MB(i/) and all remaining SBL data sets may be achieved 
if the fits are to be repeated with the low energy region 
(200< E^^ <475 MeV) excluded from the BNB-MB(i/) 
data set. 

A global analysis with BNB-MB(i/), CDHS, and at- 
mospheric constraints excluded from the fit yields a x^- 



probability of 82% and >90% compatibility among all 
experiments. 



VII. CONCLUSIONS 

We have re-examined global fits to sterile neutrino os- 
cillation models, using new data from MiniBooNE. Those 
include the final MiniBooNE neutrino mode results and 
the first, low statistics MiniBooNE antineutrino results, 
as well as first results from the off-axis NuMI beam ob- 
served in the MiniBooNE detector. 

Within a (3+1) CP- and CPT-conserving scenario, we 
have found that the data set collected by MiniBooNE us- 
ing the NuMI off-axis beam (NUMI-MB) currently pro- 
vides very weak constraints to sterile neutrino fits, due 
to large systematic uncertainties. Updated NuMI results, 
expected soon, should have a greater impact on these fits. 

Within the same oscillation framework, the Mini- 
BooNE antineutrino data set (BNB-MB(i^)) is found in 
agreement with LSND, yielding, in a combined analysis 
with LSND and KARMEN under a (3+1) oscillation hy- 
pothesis, a x^-probability of 29%, and best-fit parameters 
similar to those of LSND. Updated MiniBooNE antineu- 
trino appearance results, with almost twice the current 
statistics, are expected in the near future. 

The MiniBooNE neutrino data set (BNB-MB(i/)), al- 
though suggestive of an excess that could be described 
by a (3+1) oscillation hypothesis with a probability 
of 35%, is found incompatible with the signals from the 
MiniBooNE antineutrino and LSND results. 

The remaining null appearance and disappearance ex- 
periments (NUMLMB(i/), KARMEN, NOMAD, Bugey, 
CHOOZ, CDHS, CCFR84) and atmospheric oscihation 
data impose strong constraints to the parameter space al- 
lowed by a combined (3+1) analysis of MiniBooNE neu- 
trino and antineutrino and LSND data, excluding the 
99% CL allowed region at 99% CL. However, the con- 
straints from antineutrino disappearance experiments on 
the parameter space allowed by antineutrino appearance 
experiments (BNB-MB(i/), LSND, and KARMEN) are 
weaker. In a (3+1) oscillation framework, all antineu- 
trino experiments yield a best-fit x^-probability of 86%, 
and exclude the no-oscillations hypothesis at >5.0<j. The 
best-fit parameters are similar to those of LSND, and 
correspond to a muon antineutrino disappearance am- 
plitude of 0.35, which may be addressed by upcoming 
results from MiniBooNE and MINOS on muon antineu- 
trino disappearance. Additionally, fits to all neutrino 
experiments yield a best-fit "Probability of 47% and 
exclude the null hypothesis at >90% CL. 

Furthermore, we find that with the addition of the 
new MiniBooNE data sets, the (3+2) oscillation mod- 
els provide a much better description of all SBL data 
sets compared to (3+1) models. In the case of (3+2) fits, 
CP violation allows for a significant improvement in y^- 
probability for fits to only BNB-MB(i/) + BNB-MB(i/) + 
LSND, and fits to only appearance experiments. In the 
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case of global fits, however, the effect of CP violation is 
muted, as allowing for CP violation results in a relatively 
small improvement in the fit. The -probability for the 
best-fit (3+2) CPV hypothesis is 54%, compared to 52% 
for the CPC case. The best-fit corresponds to large but 
not maximal CP violation (^45 = I.Ttt). 

The high incompatibility among appearance and disap- 
pearance data seen in the past [1^ in the case of (3+2) 
CP-violating fits still remains with the addition of the 
new MiniBooNE results. We have shown that the in- 
compatibility is a result of the BNB-MB(i/) and CDHS 
data sets and atmospheric constraints. The compatibility 
between appearance and disappearance data with these 
three experiments excluded from the fits is significantly 
higher, at 24%. 

Neutrino and antineutrino results are also incompati- 
ble within a (3+2) CP-violating scenario, with a PC of 
less than 0.1%. The compatibility improves to 6.8% in 
the case of comparing appearance-only neutrino versus 



antineutrino results. 

Overall, allowing for mixing with multiple sterile neu- 
trino states and CP violation does not seem sufficient to 
allow incorporating all SBL experiments within a CPT- 
conserving, sterile neutrino oscillation framework. It may 
be that there is an issue with one or more of the following 
data sets: LSND, BNB-MB(i^), including the low-energy 
excess, CDHS, or atmospheric constraints; alternatively, 
theories with CPT- violating oscillations or effective CPT 
violation [ss'-'ST^ may succeed in reconciling all short- 
baseline oscillation signatures, and should be explored. 
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